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Abstract 

The spin-polarized homogeneous electron gas with densities p| and pi for electrons with 
spin 'up' (t) and spin 'down' (|), respectively, is systematically analyzed with respect to its 
lowest-order reduced densities and density matrices and their mutual relations. The three 
2-body reduced density matrices 7|-f, 7jj^, 7^ are 4-point functions for electron pairs with 
spins tt, ii, and antiparallel, respectively. From them, three functions G^^{x,y)^ G^(a;, y), 
Ga{x,y), depending on only two variables, are derived. These functions contain not only 
the pair densities according to g^^{r) = G||(0,r), ^(^(r) = G||(0,r), Qair) = Ga(0,r) 
with r = |ri — r2|, but also the 1-body reduced density matrices 71 and 7^ being 2-point 
functions according to 7^ = Psfs and fs{r) = Gss{r, 00) with s =t, i and r = |ri — r[\. 
The contraction properties of the 2-body reduced density matrices lead to three sum rules 
to be obeyed by the three key functions Gss, Ga- These contraction sum rules contain 
corresponding normalization sum rules as special cases. The momentum distributions n^{k) 
and ni{k), following from /^(r) and fi{r) by Fourier transform, are correctly normalized 
through /s(0) = 1. In addition to the non-negativity conditions ns{k), gssir), ga{r) > 



[these quantities are probabilities], it holds ns{k) < 1 and gssi^) = due to the Pauli 
principle and ga{0) < 1 due to the Coulomb repulsion. Recent parametrizations of the 
pair densities of the spin-unpolarized homogeneous electron gas in terms of 2-body wave 
functions (geminals) and corresponding occupancies are generalized (i) to the spin-polarized 
case and (ii) to the 2-body reduced density matrix giving thus its spectral resolutions. 



1. Introduction 

The homogeneous electron gas (HEG) is an important and widely used model for the 
phenomenon called electron correlation . Its model parameters are p| and pi, the 

densities of electrons with spin-up (t) and spin-down (|), respectively. Equivalent param- 
eters are the density p = Pi + Pi (from which follows the density parameter defined by 
47rr^/3 = 1/p in a.u.) and C = (Pt ~ Pi)/P' spin-polarization, is the (Wigner) radius 
of a sphere containing on average one electron. It measures simultaneously the coupling or 
interaction strength. The spin-dependent Wigner radii r^i and r^j are analogously defined 
by Airr^^/S = l/p^ and Airr^^/S = 1/pj,, respectively. This (jellium) model is relevant for the 
understanding of many effects in simple metals and semiconductors and it plays a crucial 
role in providing input quantities for approximate approaches to the many-electron problem 
of nonuniform density in solid-state theory and quantum chemistry. So the local-spin 
density approximation (LSD) of the density functional theory relies on e{rs,(), the energy 
per particle of the HEG. Corresponding LSD (and beyond-LSD) functionals E^dPh Pi] often 
treat the electron-electron interaction as a whole, without splitting it into its spin-resolved 
(11,11, and 'a' for spin-antiparallel) components. But some of the most popular functionals 
are constructed by considering spin-resolved contributions to the correlation energy. This 
makes the recent interest in the spin-polarized HEG understandable. The system under 
consideration has the two afore mentioned model parameters and (. Its ground state 
is characterized by the reduced densities n^{k), ni{k) [the momentum distributions] and 
(7^1 (r), gii{r), Qair) [the pair densities (PDs)]. Corresponding energies are t-f, [the kinetic 
energies] and f , , Va [the interaction energies] such that the total bulk energy e is given 
by e = e-f + where e-f = t| + and Cj = tj^ + with = v-^-^ + Va and f j = vn + Va- AH 
these quantities depend parametrically on and The PDs have been studied recently 
in Refs. Recent studies of n{k) for the spin-unpolarized case are in Ref. 

Q. The total (spin-weighted) g = {p]/pfg]] + {pi/ pfgn + '^{p]Pi/ P^)ga is related to 



the total momentum distribution n{k) = ^[n^{k) + ^^(A;)] (i) by the virial theorem jiol 
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The latter 



and (ii) through the 2-body reduced density matrix (2-matrix) 
means: its diagonal elements give the PDs and its contraction yields the 1-matrix, which is 
the Fourier transform of the momentum distribution. 



These mutual relations between the PDs and the momentum distributions are investi- 
gated for the spin-polarized HEG in the following. In detail its 2-matrix and quantities 
derived from it are systematically analyzed in particular in terms of the already mentioned 
contraction sum rules (SRs) and the so-called spectral resolutions. This more general 
analysis (with introductory parts in Sees. 2-6 and the main parts in Sees. 7-9) includes 
the special cases of 'no interaction' (but ( 0), 'no polarization' (but 7^ 0), and 'full 
polarization' (with C = 1 and 7^ 0) [they are summarized in Apps. 1, 2, 3, respectively]. 
Thereby four questions are considered: 

1. For the spin-unpolarized HEG it is well-known as already mentioned, that the combined 
virial and Hellmann-Feynman theorem relates the interaction energy v to the kinetic energy 
t, such that t{rs) can be calculated provided v{rs) is known from the PD g{r;rs). The 
question arises whether there are theorems for the spin-polarized HEG, such that not only 
t, but also its spin-components and can be calculated froin^ 'i>|_|^'?7jj^Ua, which are 



available because «,„ 9,, - known as function, of. and mm The answer 
in Sec. 3 is 'no' (disagreeing with a formula in Ref. |lfil|). 

2. The virial theorem is an integral relation between n{k) and g{r). Other relations arise 
from their common origin, namely the 2-body reduced density matrix. The question 'how 



these relations simplify themselves 



'or homogeneous systems ?' is answered in Sec. 8. 



3. Kimball ^| and Overhauser 3 introduced a parametrization of the HEG-PD in 
terms of 2-body wave functions (so-called geminals), which are the solutions of a 2-body 
Schrodinger equation with an appropriately chosen effective interaction (screened Coulomb 
repulsion), which describes the scattering of two electrons in the medium of all the others. 
Fn.be. developments of tb,s ,dea allows one to calcnlate tbe PD ,n excellent 

agreement with the results of quantum Monte Carlo calculations l^,!^. Because nowadays 
the PDs not only for = but also for ^ 7^ are available 0, Q,lq| , the question may be 
posed whether in this more general spin-polarized case the PDs can be parametrized again 
in terms of geminals following from a screened Coulomb repulsion. In Sec. 9 the answer 



'yes' is assumed and a corresponding analysis is presented as the necessary pre-stage of a 
numerical study. 

4. In Ref. [2^ it has been shown how the PD-normalization can be expressed in terms of 
scattering phase shifts of the geminals similar as this is known from the Friedel sum rule. 
Moreover in Ref. 2^ the Kimball- Overhauser approach is generalized with the assumption 
that the PD-geminals can be used to represent not only the PD but also the 2-matrix. 
The contraction of which gives the 1-matrix and thus the momentum distribution. Along 



this line the normalization SR of Ref. 



21| is generalized in Ref. 



221 to a contraction 



SR, which allows one in principle to self-consistently calculate the momentum distribution 
provided the afore mentioned geminals are known. For a summary cf. Ref. Q|. Again the 
question arises whether this analysis can be extended to non-vanishing spin-polarization. 
Preliminary considerations are presented in Sec. 10. 



The paper is structured as it follows. Sec. 1 defines the HEG by its Hamiltonian and 
lists the ground state quantities. In Sec. 2 the virial and Hellmann-Feynmam theorems 
are revisited. Sees. 4 and 5 present the 1-matrix 71 and the 2-matrix 72, respectively, 
with the spin-structure of 72, derived in App. 4. Sec. 6 summarizes the four PDs 
and their normalizations. Sec. 7 describes how the 1-matrix can be obtained from the 
2-matrix by contraction. Apps. 1, 2, 3 summarize the special cases of 'no interaction', 
'no spin-polarization', 'full spin-polarization', respectively. In Sec. 8 the 2-matrix is 
discussed in terms of Wick's theorem and size-extensivity, guaranteed by linked diagrams 
and cumulants. Sec. 9 presents the spectral resolutions of the 2-matrix and the cumulant 
2-matrix in terms of geminals and corresponding occupancies. Sec. 10 expresses the 
normalization and contraction SRs in terms of the scattering phase shifts which describe 
the asymptotics of the geminals. 



2. The system and its ground state properties 

The GS of the spin-polarized HEG is characterized by the average densities and pi of 
the spin-up, respectively spin-down electrons or equivalently by the average spin-summed 
density p = Pt + Pi spin-polarization C = (Pt ~ Pi)/P) thus Pt/p = (1 + 0/2 and 

Pi/P = (1 ~ 0/2- Essential parts of the phenomenon 'electron correlation' are contained in 
the two dimensionless momentum distributions n|(/c), and in four dimensionless PDs, 



namely g-^-^ir), gii{r), and g±{r) or equivalently ga = {g+ + g-)/2, Ag = {g+ - g-)/2. The 
index '±' stands for singlet /triplet, 'a' for spin-antiparallel, respectively. To get these GS 
quantities one has to solve 

H<if = Ns"^ = Ns^ (1) 

with ^ = f + y and 

Y ^ f ^fl^ J. j-^^ _ pf^^)-y^2(p2 - Pn2) :] , vi2 = — , ri2 = |ri - r2|, 
J 2! ri2 

Pi = -01-01, : P1P2 ^/'i^2V'2V^i , Pni = P5i0n(ri) (2) 

and to take the TDL Ng, Q — > 00, A^s/Q = const = ps with A^" = N-^+N^. Here 1 = (r^, (Ti), ct 
and s = t or J,, and 0i = ■0sj(r^). 'V^Kr) and ■0s (r) are the Fermion creation and annihilation 
field operators such that Ps(r) = ■0|(r)'0s(r) is the density operator for electrons with spin s 
and p(r) = pi (r) + pj^ (r) is the spin-summed density operator. Hence, Ng = J d?r ■?/'J(r)-?/'^(r) 
and N — TV^+TYj^ are the particle-number operators for spin s and spin-summed, respectively. 
The operators (: pip2 PiiPn)^ + (pn)iP2, and (pn)i(pn)2 describe the electron-electron 

( ), electron-background ( — h and -I — ), and background-background (++) interaction, 

respectively. Corresponding expectation values are the spin-density Ps(r) = (Ps(i")) for 
s =t, I, as well as the two 1-matrices ps/s(|ri — r'^^ |) = {'ipl{r[)ips{Ti)) and the four dimension- 
less PDs gss,ga,Ag defined by p^gssiru) = {■ /3s(ri)p,,(r2) :), PtPifi'a(^i2) = (pt(ri)p|(r2)), 
and p|PjA(7(ri2) = (t(^'|(ri)i/.'|(r2)i/'j,(r2'0j.(ri)). Then the spin-weighted PD g{r) is defined 
by P^g{f) — P^fl'TT(^) + plgair) + 2ptPj.fl'a(^)- The momentum distributions are given by 
ns{k) — {a\yfls)i) and n{k) — \[n^{k) -\-ni{k)\ with X]k^s(^) ~ and 2^j^n(A;) = N; the 
normalizations per particle are jfY^]^n^{k) = |(1 + C) > |, :^Z]k"'i(^) = 1(1 ~ C) < |, 
and IfEic^I^) = 1- 

3. The virial and Hellmann-Feynman theorems 

From the five mentioned reduced densities n^{k)ini{k),g^^{r),gii{r),ga{r) follow five con- 
tributions to the total energy E^T + V with T = T-f and V ^V^^ + + 2K, being 
the GS expectation values oiH^f + V with f = i\ + and F = V^T + + 2K- Now, 
the virial theorem is one relation between T and V with the consequence: if V is known, 
then T can be calculated. The question arises: Are there perhaps two virial theorems, one 



for spin | and one for spin |, which - for given PDs Qss and Qa - would possibly allow one to 
calculate the spin components T| and Tj, separately ? The answer is 'no'. To show this, one 
may define Hs = fs + Vs = Vss + Va and = + K, = Vss + K, for s =t, i, such 
that H = + Hi and E = E-^ + Ei, but there is only one virial theorem, to be derived from 
the hypervirial theorem {[E, H]) = 0, which holds for an eigenstate "if of H and an arbitrary 
operator E of the system; the commutator is denoted by [A, B] = AB — BA. With the virial 
operator F = r p and with its spin resolution rp = (rp)| + (rp)j, (rp)^ = / d?r 'iP\{y)y~^-iI)s{y) 
it is easy to show 

= -^([rp,^]) = ([P^'^]) + ([(-r|:),V']) 

Only + V ^ + V^^ depend on the system parameter VL = volume. And there is only 
one Hellmann-Feynman theorem 

The reason: There is only one MB Schrodinger equation {H — E)'$ = 0. In addition to the 
virial theorem Q one may derive the relations, cf. j^l 

= ([(rp)t, H]) = 2Tt + V11+X + m^Hi), (5) 

= ([(rp)i, H]) = 2Ti + Vu-X + m^H^), (6) 
where the expection value of 

X = U d'n£r,p,{r,)piir,)d^f. (7) 
^ J ^12 ^12 

appears. In Ref. 2q] the modified ('filtered' with the prefactor Rr/r^ = RC,/r, where R 
is the centre-of-mass and r = ri — Coulomb repulsion in Eq.((Zj) is referred to as 'spin 
polarization potential'. Note {dHg/dfl) 7^ dEs/dVt: The HF theorem does not split into 
two different theorems for spin | and spin |. So, the sum of Egs.® and (P) yields correctly 
the VT (121), but their difference gives 

d 

2AT + AV + 2X + 3Q(—AH) = (8) 

oil 



with AT = - Ti, = V|| - (K cancels), AH = - = AT + AV. 



Specification of fi: For an atom with a pointhke positive background charge, i.e. f2 — >■ 0, 
p oo, pQ — > Z (atomic hmit), all the terms with d/dVL are no longer present. Then the 
virial theorem (jHl) is simply 2T + V = Q and the other relation (jH)) is 2 AT + AV + 2X = 0, 
in Ref. j3| called 'spin virial theorem'. Contrary to this for the spin-polarized HEG in its 
TDL with A^,n ^ oo and t = T/N,v = V/N,e = E/N^nd ^Qde/dQ = r,(9e/5r„ from 
the theorems Q and Q follows the bulk virial theorem jlO| 

d 1 d d 

2t + v + r,—e = or — — r^t + — = 0, (9) 

org Ts OTs OTs 

where the derivative is related to the pressure of the system: rgde/drg = —3p/p. For the 
case 'no interaction' the virial theorem Q simply says t° ~ lAf- Also the surface virial 
theorem 2tsurf + "^^surf + {^sd/d^s + 2)esurf = can be derived [26|. But what results from 
Eq.®? One may write the sum and the difference of the last terms on the rhs's of Eqs.(jni) 
and (jBI) as 

and 

3n(AAff> = 3nAA£-3n(A)^(Aff). (ii) 

respectively. (S/Sfi)^, means differentiation with respect to \E', not H. Now, the last term 
on the rhs of Eq. fjlOj) vanishes, because of the Hellmann-Feynman theorem (^, thus giving 
the bulk virial theorem Q. Contrary to this the last term on the rhs of Eq. pij) gives an 
expression, which does not vanish. So, Eq. (fTT|) used in Eq.® gives only the trivial identity 
= 0. This is seen from 

2At + Av + rs^ = -^rlAt + —rsAv = rAj—] -— 12 



and (in the TDL) 



X 



"2 J^'^pJ 7r7 /, Y^9niR,r,0'- = (13) 



N p2 

where p^pign{R,r, ^) = (p|(ri)p|(r2)), R = (ri + r2)/2, r = (ri - ra), ^ =^(R,r), and 
goo{R, = Qaij)- So, from the HEG-PDs g^^, g^, ga as functions of r and the parameters 
TsiC can calculate only t as a function of and but not its spin components t^iti, 



separately. Thus Eq.(3.11) of Ref. |3] is wrong. But one may define a quantity y as a 
function of and ( by 

dAe 

2At + Av + r,—— = y (14) 

OTg 

(so this is not a theorem !). For 'no interaction' it is = 0, At^ ~ 1/r^, so y vanishes 
with — > 0. - In the hterature there are hints on a magnetic virial theorem if one adds to 
the Hamiltonian (0) a Zeeman term ~ 5*^ = |(iV| - iV^) Q,!^- 

4. The 1-matrices 

The 1-matrix is defined by 7i(l|l') = {'ijjy-ip]), normahzed as Tr7i = + = N . In 
terms of natural spin-orbitals V'fc(l) ^"^^ corresponding occupation numbers the spectral 
resolution of 71 is 7i(l|l') = X^fc V'fe(l)i^feV'fc(l') with < z/^ < 1 and J2k^k = N. For 
a spin-independent Hamiltonian as H of Eqs.(IH) and Q it is ipk{^) = 5s,o-i</'k(i"i) with 
(V'l(r')V'i(r)) = and (^}(r')^t(r)) = 0. Thus 

7i(l|l') = 5^^,^J(5^,,t7t(ri|r;) + (5^,,l7l(ri|r;)] (15) 

with 

j^r (16) 

The spin-traced 1-matrix is 7 = 7^ -|- 7j^ = Pifi{r) + Pif\{r). 

For a homogeneous system the natural orbitals are plane waves V5k(r) = "^^^^"^ aiid it is 
7s(r|r') = ps/s(|r — i"'!)? which defines the dimensionless 1-matrices 

fs{r) = jrY.^s{ky''\ (17) 
k 

displayed in Fig. 1. The Fourier transforms of /s(?") give the momentum distributions ns{k) 
with < ns{k) < 1 and with a large and a small Fermi sphere: {k-p^Y = Q'^'^Ps, ^ft > ^fj,- 
With kl = SvrV it is kp^/kp = (1 + C)^/^ and kpi/kp = (1 - C)^/^- /.(O) = 1 makes the 
ns{k) correctly normalized: X]k^.5(^) ~ The deviations of ns{k) from and 1 (its 

non-idempotency) measure the correlation strength: X]k[^«(^)]^ ^ ^s- The spin-summed 
momentum distribution is n{k) = ^[n-^{k) +ni{k)]. Whether the asymptotics of the ns{k) is 
spin-independent has to be checked. If so, it should hold ns{k 00) = (2a;o/37r)^5fa(0)//c^, 
ao = (4/97r)V3r, Q,Q- 



For the special cases 'no interaction', 'no polarization', and 'full polarization' cf. Apps. 
1, 2, and 3, respectively. 



The recently available momentum distribution n{k, Vg) of the spin-unpolarized HEG {k 
measured in units of kp) can be used to approximately construct the momentum distri- 
butions of the spin-polarized HEG for small values of (: 

^T^^'"-^)^^((rT7F^'(rTtF^)' ^^^^'"-^^^^ ( (1-0^/3 ' (1-0^3 )- (^^^ 

This allows one to calculate the above mentioned spin components t^{rs,() S'lid ti{rs,C) of 
the total kinetic energy t{rs,() = t^{rs,C) + ^j,(^s;C) ^-^id to control it by comparing with 
t{rs, 0, calculated from the PDs with the help of the virial theorem Q. Results for ( = 1/3 
and Ts < 10 are displayed in Figs. 2 and 3. They show that the approximation (|T8|) affords 
the correct t{rs,C)- 

Related quantities derived from /|, fi or n^, are the 1-matrices (referred to as non- 
idempotency matrices) 

Psirw) = I dh2pfs{n2)fs{r2v) = ^ X]K(A;)]2e^'^^"' , /3.(0) = 6., 

/?a(rnO = ^ / d\2Ph{r,,)f^{r^y) = l^Y.n^{k)n^{k)e^^^^^' , /3,(0) = 6,. (19) 

The product n-^{k)ni{k) makes the k-summation to run over three different regions: The 
inner sphere A; < (1 — CY^^i the shell (1 — C)^^'^ < A; < (1 + C)^^'^) ^^nd the outer region 
> (1 + Cy^^- The quantities 6^ and ba are related to spin-generalized Lowdin parameters 



k 



Ca = l-^K = l--^5^nt(/^)ni(A;), (20) 

which vanish for 'no interaction', cf. App. 1. For 'no polarization' (cf. App. 2) it 
becomes n^{k) = ni{k) = n{k) and therefore hg = ha = h, Cg = Ca = c, c = 1 — 2b, and 
b = jf Ylkl'^i^W ■ ft already P.-O. Lowdin who has asked for the meaning of Tr7^. 
The answer is nowadays: 1 — Tr7f/A^ is the normalization of the cumulant 2- matrix and 
one possible index of the correlation strength. As seen from c = jf^]in{k)[l — n{k)], this 
expression is particle-hole symmetric 31[. In addition to the above defined non-idempotency 



matrices, the 1-matrices a(r), a{r) appear, when contracting the 2-matrix, cf. Eq.()36p. 
They cause (for spin-antiparallel electron pairs) the difference between singlet and triplet 
PDs and describe fluctuations, which mimick spin-flip processes. This means the following: 
One may consider a spatial part (region, domain, or fragment) u of the system, say a sphere, 
which contains on average up electrons. One may furthermore ask for the probability of 
finding A^^^ 7^ t^P electrons in u. This can be realized by a certain number N^-^ of spin-up 
electrons and a certain number N^^ of spin-down electrons. These numbers can change 
under the constraint N^^^ + N^^^ = N^, what looks effectively as if electrons in uj had 
made spin-flips, although such fluctuations are based only on electrons (with spin-up and 
spin-down) leaving or entering the region uj. In general, for particle number fluctuations in 
spatial regions cf. Ref. j^^. - For 'no interaction' it is a^{r) = -j^f^ir) and a°(r) = 0, cf. 
App. 1. 



5. The 2-matrices 

The 2-matrix 72(1|1', 2|2') = {■ijj\,ipl,-ilj2ipi) is with N^^-ipi = ipi^N^^ — ^o-io-a), hence Nipi = 
ipi{N — 1), normalized as Tr72 = A^(A^ — 1) with Tr - ■ ■ = J dld2- ■ ■ |i/=i 2'=2- For given 
positions r^, r'^^, r'g there are 16 spin matrix elements. But only six of them are non- 
zero, the other 10 matrix elements vanish. With the short-hand notation {a[a2cr2<^i) = 
(-0^, (r'^)V'^/ (r2)V'o-2(r2)^o-i(ri)) the non- vanishing elements are 



7tt = 


(TTTT), 


Tr7|| = 




111 = 


( j, J, j, j,) 5 


Tr7ii = 




= 


(Tiit), 


Tr7^ = 




ll = 


(itti), 


Tnl = 




ll = 


(itit), 


Tr7f = 


-Na, 


ll = 


(Titi), 


Tr7^ = 


-Na, 



whereas e.g. (ttii)) (TTTi) vanish. Each of these matrix elements depends on ri|r'j^, r2|r2. 
The first two lines concern spin-parallel electron pairs and the next four lines describe 
spin-antiparallel electron pairs. The properties of these matrix elements 1)211) are listed in 
Table I together with other related matrix elements for the spin-antiparallel electrons and 
with the corresponding spin-matrices defined below in Eq. ()22|) . The operators therein are 
defined as it follows: D creates the diagonal elements with r'^^ = ri and r2 = r2, P makes 



the permutation ri ^ r2, Q effects the hermitian conjugation with r'^^ ^ ri, r'g ^ r2, and 
c.c. (= complex conjugation), and C causes the contraction with ^ V2 and j d^r2, 
where 1-matrices appear, namely 7f = P|/t, 7j, = Pj,/j,, and pas, pa4. For the position 
matrices 7, Tr means r'^ = ri, r'g = r2, and J (Prid^r2. For the spin-matrices 5, Tr means 
o-'i = ai, 0-2 = 0-2, and Z1<ti,2- ^^^^ 7a = (7- - 7+)/2 and = {g+ - g^) /2. Columns for D 
making r'^^ = r2, Y'2 = ri and similarly for P, Q, C are not shown. - In the column D7, i.e. 
D7(ri|r']^, r2|r2) = 7(ri|ri, r2|r2) the PDs g^-^, g^, ga, and Ag are defined. They depend on 
^12 = ki — r2|, whereas all the position matrices depend on ri|r'^, r2|r2. Note = D7^, 
which holds because of (V^|(ri)'?/;|(r2)-?/'j(r2)V^l(ri)) = (■?/'|(r2)-?/'|(ri)-?/'|(ri)-?/'|(r2)) and 
r2i = ri2. It is similarly D7^ = D7^, but for the off-diagonal elements it is •yl 7^ 7^ and 
7a 7^ 7a- Note further, that the spin-antiparallel 2-matrices 7^, 7^, and 7^, 7^ are indeed 
hermitian, but have not the Pauli-principle property P7 = ±7. 

To calculate V^s and Va we need only the diagonal elements D'jss and -07^^. Also for the 
contractions only C'~fss and C7^'^ are needed. So the question arises: what is the physical 
meaning of the spin-exchange matrices 7^''^, which information they contain ? Because of 
D7^ = — D7^, D7^ = — D7^ they contain in principle the same PD-information as 7^'^. 
The 2-matrices 7^'"^ mimick spin-flip processes such that the total spin is conserved. Thus 
they play a role only if the particle interaction would be spin-dependent, e.g. different for 
singlet and triplet spin-states. The normalization of 7^'^, Tr7^'^ = —Na{rs,C)j does not 
appear in the normalization of the total 2-matrix 72, because of Tr5^'^ = as shown be- 
low in Eq. ()22|l . For the special cases 'no interaction' and 'no polarization' cf. Apps. 1 and 2. 

With the spin-matrices (corresponding to the six matrix elements jss, la, ' ' , iti namely) 



5^{ai\a[, CT2 


a'2) 


= 5<7l 






-'2 




Tr(5-|-| = 


= 1, 




4) 


= 5<71 




^02 


-'2 


^<TlA^02A^ 


Tr5u = 


= 1, 




a'2) 


= 5<71 




^a2 


-'2 




Tr5^ = 


1, 




a'2) 


= 5<7l 




^a2 


-'2 


^tTl,l^(T2,h 


Tr5^ = 


1, 




^2) 


= 5<71 




^CT2 


< 




Tr5f = 


0, 


'^l(^lki>^2 


a'2) 


= ^ai 




^a2 






^< = 






(note D(5i^'^ = 0) the total 2-matrix 72 is given by 

72 = ^ttttt + ^laii + ^bl + ^hl + ^hl + ^hl 

Tr72 = A^t(A^t - 1) + Ni{Ni - 1) + N^N^ + N^N^ + + = A^(A^ - 1). (23) 
From Eq.(j2SI) also follows the spin-traced 2-matrix 7|| + 7^ + 27^, where 7a = |(7a + 7a)- 

In view of the spectral resolution of 72 (respectively of its position matrices) in terms 
of antisymmetric spin geminals, the drawback of Eq. lj^ is that 7a5 ■ " " ^it eigen- 
functions of the permutation operator P as seen from Table I, column P7. Only the linear 
combinations 7± have the eigenvalues ±. So it is useful to rewrite Eq.(j2ni) by appropri- 
ate (the Pauli symmetry achieving) linear combinations. In the first (intermediate) step 
with the definitions of Table I for 5a-, 7a etc. (the last but one block of lines) it results the 
parallel/antiparallel representation 

72 = 5||7|| + 5^7^ + 5a7a + ^'aia + ^a7a + Kl'a^ 

Tr72 = A^t(A^t - 1) + Ni{Ni - 1) + 2N^Ni + + + = N{N - 1). (24) 
In the next step using 5±,7± etc. of Table I (the last block of lines) it turns out 

72 = [5tt7tT + hllll + '^+7-] + ^-7+ + [^-1+ + KlA^ 
Tr72 = [N^{N^ - 1) + N^{N^ - 1) + {N^N^ - Na)] + {N^N-^ + Na) + [0 + 0] = N{N - 1). 

(25) 

This is the triplet /singlet representation of the 2-matrix 72, because the first three terms de- 
scribe the spin-triplet, the next term the spin-singlet, and the last two terms a singlet /triplet 
mixing, cf. App. 4. Note, that in the trace the terms ±A^a and in the diagonal elements 
the terms ±piPiAg cancel each other. The symmetrized, respectively antisymmetrized spin- 
matrices used in Eq. ()25j] are 

5±(0-l|0''l,0"2|0'2) = 2^^ai,<T[ Sa2,ai,±^au<T!2 ^a2,a[){Sai^Sa2l + Sa^l6a^^), TtS± = 1, 

6±{ai\a[,a2\a'.2) = ^{^ctuct[ ^^2,^'^ T S^.y^ S„^^„'J{6„,^6„.,i - 6^^i6„^^), Tt6± = (26) 
and the corresponding position matrices are 

7± = lill + 7a T 7a T 7a), Tr7± = N^N^ ± Na, 

7± = ^(7a - 7a ± 7a T 7a), Tr7± = 0, (27) 



as it follows from Table I. 

Hermiticity: Whereas j^-^, 7^, and 7-1- are hermitian, the last two terms of Eq.(j2Hl) are 
hermitian only summed together because of Q6± = 6^ and Q'y± =7^- - As above the 
spin-traced 2-matrix is again 7|| + 'y^ + 27^, because of 7+ + 7_ = 7a + 7a = 27a. - Note 
that the interaction energy V consists of a triplet contribution V^^ + Vn + V- arising from 
7ss, 7_ and of a singlet contribution stemming from 7_|_. 



6. The pair densities 

For the diagonal elements 72(1|1,2|2) = P2(l,2) = p'^g{l,2), from Eq.^ and Eq.(j2SI) the 
two equivalent expressions (r = ri2) 



or 



'^aiT<^a2TPT^TT('") + ^a^iSa2lpl9ll{r) + lK,-a2plPl9-{r) 



+ ^^a,,^a2plPl9+{r) 

(28) 



result. The two dimensionless spin-parallel PDs giiir), g^i^r) are (with r'^^ = ri, r2 = r2) 
defined and with Tv'-fss = Ns{Ns — 1) normalized by 

Pti/nW =7TT(ri|ri,r2|r2), j d^r p[l ~ g^^{r)] = g^^{oo) = 1, 

Pl9ll{r) = 7u(ri|ri,r2|r2), j d^r p[l - g^{r)] = ^u(oo) = 1. (29) 

The spin-antiparallel PDs ga{f), g±{r) are defined and with Tr7a = N^N^, TT'-f± = N-^N^±Na 
normalized as 

PTPi^a(^) = 7a(ri|ri,r2|r2), j rfV p[l - ^a(r)] = 0, ^a(oo) = 1, 

f 4a 
PTPi^±('") = 7±(ri|ri,r2|r2), y rfVp[l-^±(r)] =^Y— ^±(00) = 1, (30) 

where 4a/(l — C^) equals 2/{l + () for 'no interaction' and 2 for 'no polarization', cf. Apps. 
1 and 2. (7+ and g^ are the singlet and triplet spin-antiparallel PDs, respectively. (Here 



they are normalized as g±{oo) = 1, whereas in Refs. [2l|, |2^ it is g±{c)o) = 1/2.) For 
spin-independent particle interaction the PDs (?+ and g_ are equally weighted, so only their 
sum g+ + g- = 2ga appears (and not \{g+ — (?-) = ^g, which is the singlet-triplet splitting 
of the PDs g^ and g^). Spin summation yields the spin-weighted PD g{r) = J2a 5'(1)2), 



which is also given only by Qss and ga- 

9ir) = {^?g,,{r) + {P^fg,,{r) + 2P^g,{r) or 
1 - 9{r) = {^f[l - gu{r)] + (^)^[1 - gn{r)] + 2^[1 - g^{r)] (31) 

r r r 

with pt/p = (1 + 0/2 and pjp = (1 - 0/2 and / <Pr p[l - g{r)] = 1, (7(00) = 1. 
'On-top', i.e. for r = 0, the coalescing cusp theorems 



hold 



HQ- 



3 / 4 

^^^(0) = 2"o^"s(0)' ^l(O) = ao^a(O), «o=(^l (32) 



7. The contraction SRs 

The contraction of the spin-parallel 2-matrices, 

j rfV27s.(ri|r;,r2|r2) = 7,(rJr;)(iV, - 1) or 

j d^r2[psls{ri\r\) - 7ss(ri|ri, ralrg)] = 7s(ri|r'i), (33) 

makes the spin-parallel 3-point functions 7ss(r]^|r'^, rgjrg) to appear. Because of translational 
and rotational invariance (six conditions) this function depends only on the three distances 
fiv = |ri — r'll, T12 = |ri — r2|, and ri'2 = \r[ — r2\, besides the direction of ri2 can be 
averaged: 

Gssirir,ri2) = J ^^^^^^^(rii/, ri2, \rir - r^l), PsG'ss(^ii', ^2, ^'2) = 7,s(ri|r'i, r2|r2). 

(34) 

Gii{rii',ri2) and G'n(rii/,ri2) are special 2-matrices yielding with r[ = ri the spin-parallel 
PDs: gss{ri2) = Gss{0,ri2). The contraction SRs in terms of these 2-matrices 

G'ss(rii/, ri2) are rewritten as 

d^ri2 Ps[fs{riv) - Gss{rir,ri2)] = /s(rii'), G^sirir, 00) = fsinr). (35) 

So the functions Gss contain not only gss{r), but also /s(r) and with rn/ = the SRs (jHKjl 
yield the normalization (pHjl of gss{i^)- 



The contraction of the spin-antiparallel 2-matrices 7-1- yields 

j rfV27±(ri|r;, rzira) = ^[N^-f ^{ri\r[) + N^-fi{ri\r[)] ± pa{ru') or, 
j d^r2 |^h7T +PT7i](ri|r;) -7a(ri|r;,r2|r2)| = ^pa(rn/). (36) 
Corresponding 2-point functions are defined by 

With this definition of G± the contraction SRs can be rewritten as 

I d'rup[f{nr)-G±irn',n,)] = T^^^, G±(rn', 00) = /(mO (38) 

with /(r) = |[/t(r) + /j,(r)]. The functions G± contain the spin-antiparallel PDs g±{ri2) = 
G±{0,ri2). Thus the SRs ()38p contain with rn/ = the normalizations ()30|) of g±{r). Note 
that the 1-matrix a{rii/) determines the normalization of the PDs g± with a(0) = a. - For 
Ga = {G+ + G'_)/2, which equally weights the singlet /triplet contributions, it holds 

j d\i2 p[f{riv) - Ga{riv.n2)] = 0, Ga{ru', 00) = /(mO- (39) 

In summary, the two functions Ggsirw ,ri2) not only contain the two spin-parallel PDs 
gss{r) = Gss{0,r), but also the two 1-matrices /s(r) = Gss{r, 00). And the function 
Ga{rii' ,ri2) contains the spin-antiparallel PD ga{r) = Ga{0,r) and besides it controls 
/ = ^[/| + fi] with |[/t(r) + fi{r)] = Ga(r, 00). A more refined description distinguishes 
between the singlet function and the triplet function G_. So there are six SRs, namely 
the four contraction SRs (jH^j) and (jHHj) and the two momentum-distribution normalizations 
fs{0) = 1. 



8. Wick's theorem, linked diagrams, and the cumulant expansion 

If the Coulomb repulsion is treated as perturbation (with many-body S'-matrix theory and 
Wick's theorem), then the 1-matrices 7| and 7j are given by linked Feynman diagrams only, 
as sketched in fEq.(fT3j) shows with 71 ~ 5o-i,o-2- there is no spin-flip along an electron line) 



+ 



+ 



(40) 



Contrary to this the 2-matrices are given by both unhnked (or disconnected) and hnked 
(or connected) diagrams. If all the unlinked diagrams are summed up, products of the 
1-matrices appear, e.g. 

-fss = {ssss) = 7s(ri|r'i)7^(r2|r2) - Js{ri\r2)js{r2\r[) + {ssss)c, (41) 

where {ssss)c = —Xss is given by the sum of all the connected diagrams of 



Xss = (^l(r;)V^J(r2)V^s(r2)4(ri))f 



+ 



+ 



(42) 



The index 'c' means 'connected' or 'cumulant'. It is similarly 

7a = (Tiit) = 7T(ri|r;)7i(r2|r'2) + (TiiDc 

7a' = (itit) = 7T(ri|r'2)7i(r2|r;) + (itiDc (43) 
where (TiiT)c = — Xa ^'^^ (iTiT)c = — Xa mean the linked diagrams of 



Xl = (V^j(r;)V^|(r'2)^i(r2)^T(ri))c = ! + ■ ■ • , 

» ^ » ^ » 



Xl = (^j(r;)V'}(r'2)V'^(r2)V'T(ri))c = ! + " " • • (44) 

, ^ I ^ , 

The products of the 1-matrices 7s = psfs may be called generalized Hartree-Fock (HF) 
parts, where the term 'generalized' means 'built up from the true, hence non-idempotent 
1-matrix /s(r)'. The HF parts of the six position matrices ()2H1 together with eight other 
related 2-matrices and their properties are listed in Table II. The matrices Xss, Xa; xl etc. 
are called cumulant 2-matrices, they are listed in Table III. Eqs. fl41|) and (pSj) are the first 
steps of the more general cumulant expansion js^ . Per construction, the cumulant matrices 
are given by linked diagrams, what makes them size-extensively normalizable. Contrary 



to this the normahzations of the 7-matrices and the 7^^-matrices contain terms ~ A^^, cf. 
the columns 'Tr' of Table III with those of Tables I and II. In the column C7^^ of Table 
II the already above introduced 1-matrices Psirw) and Pairw) appear, as it is obvious 
from Eq. ()19j) . From the diagonal elements D7^^ follow the PDs g^^{r) = 1 — |/s(^)P with 
/ d'r[l - gffir)] = f 6, and g^^ir) = 1 as well as A^?«^(r) = Mr)f^ir). 

The diagonal elements of Xss and xl Xa S^^^ ^^e cumulant contributions to the PDs: 
P%s = Dx., / d'r phUr) = ^(1 - ^bs), gssir) = 1 - |/,(r)|2 - h,{r), 



P]Plh± = D^i'^ ph±{r) = Tia-ba),g±ir) = iTl^ir) - Pa{r)]- h±{r), 
p^PlK = Dxa, j d^r pha{r) = 0, ^^(r) = 1 - - hair). (45) 
The contraction of the cumulant 2-matrix: With the definitions 

plHss{rir,ri2) = J J d^r2Xss{ri\r\,r2\r2), 

j c?^'^2X±(ri|ri,r2|r2), (46) 

and Ha = + H_)/2 the contraction SRs are 

I d^r 12 pHss{r r 12) = — - — /3s(rii')], 

J Ps Ps 

f p^ 

/ d^ri2 pH±(rii>,ri2) = T [(^{rn>) - /3a{ru')], 

J P1P1 

j d'ri2pHainr,n2) = 0. (47) 



The advantage of the if-functions: in perturbation theory they are given by linked diagrams, 
which are size-extensive. They contain the cumulant PDs hss{r) = Hss{0,r), h±{r) = 
H±{0,r), hair) = Ha{0,r). To get also the 1-matrices fg (or equivalently the momentum 
distributions rig) one has to write the first equation of the SRs pUj) as 

j d'r,2pHUru',ri2)= ^ ^ n,(A;)[l - n,(fc)]e^'^'-"' (48) 

and to Fourier analyse the Ihs according to 

j d^u pHUr,r,r,2) = (^j-^ ^ m,(A;)e*^^"' (49) 



defining functions ms{k). Tlien one lias to solve the quadratic equation 
n^(fc) — ns{k) + ms{k) = 0. As already mentioned, one way to get the PDs and the 
momentum distributions is to perturbatively calculate linked diagrams for Hgs and H±. 
Another way is to solve an effective 2-body Schrodinger equation for the geminals, 
which parametrize the 2-matrices 7.0^ and 7-1- like this has been done successfully for the 

spin-unpolarized HEG G, Q, Q, Q- 

This way is sketched in the following for the 

spin-polarized HEG. 



9. The spectral resolutions 

The four 2-matrices 7^5 and 7-1- are needed for the key functions Gss and G±. They have 
the correct Pauli symmetry and are hermitian. So they can be diagonalized in terms of 
symmetric/antisymmetric geminals and corresponding occupancies. This is similar to the 
well-known natural representation of the 1-matrix. For the HE parts this diagonalization 
can be done easily with the help of the spectral resolution or natural-orbital representation 
(fTTjl of the 1-matrices /s(r) in terms of plane waves (because the system is homogeneous). 
Thereby the geminals 

^e'^^^e''^'- with R=^(ri + r2), r = ri-r2 (50) 

appear together with the weights 

fis(K,k) = ns{ki)ns{k2) , /ia(K, k) = n|(A;i)n|(/i;2) 

with ki = + k , ks = - k. (51) 

R is the centre-of-mass coordinate and "^e'^^ describes the centre-of-mass motion (which 
is a free-particle motion, because the system is homogeneous), whereas r is the relative 
coordinate and (/^'^(r, k) = -^Q^^'^ describes the relative motion of the HE parts. Starting 
with the expressions for 7^^ and 7^?^ of Table 11, the following spectral resolutions of the 
HE parts result 

7?/(R|R',r|rO = 5^/x.(K,k)v.° (R,r)v,°*(R',r'), 

K,k 

7r(R|R',r|r') = 5^/.4K,k)^° (R,r)^°*(R',r') (52) 

K,k 

with 95° (R, r) = -^e^^^-^ i'jL(A;r)F£(efe), where the dependence of K, k is not exphc- 
itly denoted and the angular momentum expansion of a plane wave is used. Note L = {I, mi) 



and jL{kr) = ji{kr)YL{er). So the spectral resolution ()52|) can be written in terms of free- 
electron geminals ji{kr) and occupancy matrices /i|^^,(K, /c) and fi1]^,(K,k). The normal- 
izations 

ki,2 k 

TtjT = E'"«(K>k)[l±5k„kJ=iV^iV^±E^T(fcKW=A^TA^i±A^&a (53) 

ki,2 k 

follow from J ^e''^'" = (5k,o in agreement with Table II. 

From the spectral resolution ()52|1 of the HF parts one may derive the spectral resolution 
of the total 2- matrix by replacing the the free-electron geminals or ji{kr) by interacting- 
electron geminals !f± or Ri{r, k): 

7,,(R|R', r|rO = ^ /^«(K, k)(^_(R, r)^*_(R', r'), 

K,k 

7±(R|R',r|r') = E/^-(K,k)¥^±(R,r)y.i(R,r'). (54) 

K,k 

(The general structure (jMj) in terms of eigenfunctions and eigenvalues, which diagonalize 
and 7-1- , is correct, but the use of Eq. ()51|) for the geminal weights, which is right in Eq. (j3^ . 
may be here in Eq.(jHH) only an approximation.) With the angular-momentum expansion of 
(/9±(R, r) the PDs Qss and g± are given in terms of the afore mentioned geminals -R/(r, k) and 
occupancies /i|^^,(K, fc) and /i2L/(K, fc). Their diagonalization causes a (K, A;)-dependent 
L-mixing with L-mixed geminals R^{r,k) and occupancies /x^. Calculating the PDs with 
the non-spherically symmetric functions |/2^(r, fc)p, their non-sphericity is averaged out 
by summing over A with /i^. The geminals Ri{r,k) are the scattering state solutions of 

2 

a 2-body Schrodinger equation with an effective repulsion v±{r) = ^ + which may 
be in general a non-local interaction potential, but a local one possibly can be a good 
approximation. Whithin a PDFT the ellipsis result from a functional derivative of the 
kinetic energy as a functional the PD js^. 

Eq.()54|) means for the cumulant matrices 

X..(R|R',r|r') = E/^^(K,k)b°(R,r)^°*(R',rO - V5-(R,r)(^*(R',r')], 

K,k 

X±(R|R',r|r') = E/"a(K,k)[y.°(R,r)^°*(R',r') -(/.±(R,r)¥;i(R',r')]. (55) 

K,k 



Because the cumulant matrices Xss and x± ^ire hermitian, they can be represented 
alternatively in their 'own' spectral resolution with cumulant weights /i^ ^ and cumulant 
geminals Ri being the bound state solutions of another 2-body Schrodinger equation with 
an effective attraction f^(r), possibly f^(r) = —v±{r). 



10. Scattering phase shifts 

The afore mentioned geminals Ri{r, k) are scattering states with an asymptotical behavior 

Riir,k) ^ ^sm[kr-l^+r]i{k)] (56) 

defining scattering phase shifts r]i{k) [v±{r) is assumed to decay stronger than 1/r^]. For 
the case 'no polarization' the PD normalizations 

j p[l - g.^{r)] = 2, y d'r p[l - g,{r)] = (57) 



can be equivalently reformulated as the Friedel like SR 

± r°° 1 

Y,^^ dk[-fi'{k)Uk) = ±c, p(A;) = -5^n(A;i)n(fc2) (58) 



2 



K 

with ki^2 = ||K±k|. Generalizing this to the case of non- vanishing polarization, one should 
expect e.g. for the spin-parallel PDs 

-V / dk[-fi'^{k)]r]i{k) = -c„ fis{k) = —y2n,{ki)n,{k2) (59) 

as an equivalent reformulation of the normalization (jSHI)- 3k'^fis{k) is the probability of 
finding two electrons with spin s and momenta ki and k2 with the half momenta difference 
/c = ||ki — k2| js^. Similar SRs should hold for electron pairs with antiparallel spins. Also 
the contraction SRs (29) of Ref. [2^ can be generalized to the case of non-vanishing spin 
polarization. 



Summary and outlook 

Based on the spin-structure of the 2- matrix [Eq. ()25p ]. its contraction properties ()33|) and 
flHUj) lead to three two-variable functions Gii{x,y), G^(x,?/), and Ga{x,y), which are key 
functions, because they contain not only the PDs according to gssi^) = G'ss(0,r) and 
ga{r) = Ga{0,r) with gss{0) = 0, ga{0) < 1 and gss{r), ga{r) > 0, but also the 1-matrices 



according to = Gss{r, oo) with /s(0) = 1. These key functions Ggs and Ga have to obey 
the contraction SR s (1351) and (jHIIj) . They can be approximately calculated (i) perturbatively, 
cf. Sec. SandRef. ffo. (i^ bpl.e»luiio„ of an effective 2-body e<,uat.on according to the 



SRs of Refs. 



21 



Kimball-Overhauser approach |l8l . ll9L l2r| . The Friedel like normalization and contraction 



can be generalized to the case of non-vanishing spin-polarization. The 
general analysis presented here has to be completed by numerical studies. Also the question 
should be studied, whether conclusions can be drawn for the afore mentioned key functions 
and geminals from the hierarchy of contracted Schrodinger equations |13l, the Kummer 

n n 

variety [l^jhnear inequalities for density matrices [la], the P/D, Q, and G positivity 
conditions |4^. Future work should concern also the generalization of the fluctua- 
tion analysis of the spin-unpolarized HEG to the case of non- vanishing spin-polarization. 
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Appendix 1: The case of 'no interaction' 

Because of T ~ and V ~ this case is realized for — ^ 0, what effects the following 
simplifications, where N^/N = (1 + 0/2, NJN = (1 - C)/2 and kp^/kp = (1 + C)^^^ 
kpi/kp = {1 — (Y^^, and wavelengths are measured in units of kp, lengths in units of kp^, 
and energies in a.u.: 

n?(A;) = 0((l + C)^/^-fc), n'^ik) = 9iil-CY/'-k), (60) 

There are two Fermi spheres, a large one for spin | and a smaller one for spin |. There are 
three regions of different occupancies: ti in the inner sphere, t in the shell, and outside the 



shell. Because of the idempotency = n^{k) and the projection n'^{k)n^{k) = n°(/c) 

from Eq. ()19p it follows 

= = &° = (63) 

All the quantities and relations of Table II hold [with the replacements HF — > and /s(r) — »• 
fg{r)] also for the ideal Fermi gas including Eq. ljU^ and 

44(0 = «°(r) = ^^/f (r), a°(r) = 0. (64) 

It follows a° = (1 — C)/2, what can be proved also in the following way. Starting with the 
definition a = —Tvy^'^^/N one has to evaluate Tr7^ or Tr7^ for 'no interaction', which allows 
with Eq. (p]|) to calculate the traces in the k space: 

ki,2 ki,2 k 

Thereby a|k2^jk2) should be a state with a hole at k2 in the |-sphere (hence k2 < kp^) and 
an additional particle at k2 in the j-sphere. This is not possible: ^ojkj = 0. From Tr7^ 
follows the same result: 

Tr7a = Zmki«lk2«Tk2«iki) = Eki,2(«Tki«iki«ik2«Tk2) - Ek(«TkV)- (66) 

ki,2 

In the first term, with a^^^ ^ i^ole in the f-sphere can be made only in the shell 
kpi < k2 < kpi between the two Fermi surfaces. Because k2 is thus outside the |-sphere, 
the I places in the shell are empty and can be occupied: a^ii^^a^^^^) ^ 0. With ki = k2 and 
Efc,,<fc<fc,, 1 = iV^ - iV^, again a" = N^jN = (1 - C)/2 turns out. 

The PDs for spin-parallel pairs are 



9\i{r) = \-\fl{r)\\ J p[l - gl,ir)] = gl,ioo) = I. (67) 
The PDs for the spin-antiparallel pairs are 

^7°(r) = 1, / d^r p[l - g'^ir)] = 0, g',{oo) = 1, 

g'a'ir) = -Re/f (r)/°(r), / d'r pg'^{r) = ^, g'^^oo) = 0, 

^0 (r) = 1 T Re/f (r)/°(r), / d'r p[l - gl{r)] = g^{oo) = 1. (68) 



Finally the spin-weighted PD is 

l-Ar) = (^)'l/TV)r+(^^)'|/fWP (69) 
with / d^r p[l — g^{r)] = 1, g^{oo) = 1, besides Ag^{r) = g'^{r). 

Appendix 2: The case of 'no spin-polarization' 

The simplifications for 'no polarization' are: C = 0, consequently = Ni = N/2, and also 

PT = Pi = P/2, 7t = 7i = 7/2, /t = /i = /, 7 = P/, «T 
position matrices are 

5p = + = (Jaiai'^aa^^aiaa, (TTTT) = (iiii) or 

Sa = 5l + 5l = (Tiit) = (itti) or 

5i = 5a + 5a = 5.x<.^5<x..i5.i,-.., ant) = (Tin) or 

An immediate consequence is 7a = 7^ = 7± = 0. Besides from Table I it follows 7^^ + 71^^ = 
7^^ and the Feynman diagrams (42) and (44) lead correspondingly to Xa + x'a — Xpy because 
the spin-directions t and | are equally weighted. (Note also, that 7^ and 7- = 7a + 7a 
have the same symmetry property, namely changing sign under the replacements ri ^ r2 
or r'^ ^ r2.) So it is 

7a + 7a = 7p or 7;^ = 7p - 7^, or 7_ = 7^ (71) 

(for Vs = this can be proved directly). Eq. (ffT|l says that for C = the 2-matrix 72 consists 
of only two components, namely 7^ and 7a or 7+ and 7_. Indeed the parallel/antiparallel 
representation ()24|) takes the simple form 

Tr72 = 2.^(|-l)+2.(|y = A(A-l) (72) 

and the equivalent triplet /singlet representation simplifies as 

72 = {6p + 5+)7_ + 5_7+ 
Tr-. = 34(f-l)+l4(f + l)=iV(iV-l). (73) 



ni = n. In this case the spin and 



N (N 

7tt = 7u = 7p, Tr7p = y I y " 1 

7a ~ 7a = 7a5 Tr7a = ^y^ ' 

7a^ = 7a'^7l, Tr7; = -y. (70) 



Its equivalence with Eq. (f7^ follows from Eq. (f7T|) and 7^ = (7_|_ + 7_)/2 and 6± = {6a±S'^)/2 
of Table I. The spin structure of Eq. ()73|) agrees with the result of the RPA-like treatment 
for small [Is^. Note Tr7-|- = ± !)• The spin-matrices appearing in Eqs. lff^ . (ffHjl are 

(In the last line it holds |(---) = |(' • ■ )^cri,-o-25 because of 1 = ^o-i.o-a + ^0-1,-0-2 and 
■ ■ )<^o-i,o-2 = 0.) The arrows indicate the diagonal elements, which give the spin-dependent 
PD D72 = P2 = P^g and the spin-summed PD g{r) = ^ g{l, 2) equivalently as 

0"1,2 

ai^, 2) = l[Sau^29p{ri2) + Sa,,-a2ga{ri2)], g{r) = ^[gp{r) + ga{r)] or 

^(1,2)= l[l{l + 6^,,,,)g_{n2)+'^S^,,-.29+iri2)], 9{r) = \[3g.{r) + 9+ir)], (75) 

in agreement with Eq. ()31|) . From the above relations between the p, a and the it components 
of the 2-matrix follow corresponding relations for the PD: 

7p = 7_, 7 = ^(7+ + 7_) or 7+ = 27„ - 7p, 7. = 7^, 

9p = 9-,9a = ^{9+ + 9-) or g+ = 2ga - gp, g^ = gp. (76) 

Whereas in general Ag = {g+ — gJ)/2 (= singlet-triplet PD splitting), for C = it is 
/S.g = ga — 9p (= antiparallel-parallel PD difference) or gp = g — Ag/2, ga = g + Ag/2. So, Ag 
may be addressed as magnetic PD. The contraction of Eq.((7T)) yields f /(y — 1) = ^f^ — pot 
or a{r) = /(r)/2, thus for r = it follows a = 1/2. In the contraction and normalization 
of the HE and cumulant parts appear Ps = Pa = P, thus hg = ha = h and Cg = Ca = c, 
h = (1 — 6)/2. The case C = is summarized in Table IV. 

The PDs are defined and normalized according to (with r'^^ = ri, r'g = r2, r = ri2) 

(0%±(r) = 7±(ri|ri,r2|r2) , j dh p[l - g^{r)] = t'^, g ^{00) = I. (77) 

With the functions 

(0'G±(rn^ri2) = j '^^^^{v,\r\M^2) (78) 



^(Ti,(T2 ; 



+ 



■'(Tl, — (72 • 



(74) 



the contraction SRs are 

j d\i2 p[/(riiO - Giim', ri2)] = T2f{nr) , G±(rir, oo) = /(mO- (79) 
From the spectral resolution ()54j) it follows 

Gi{nr,n,) = I ^4E^(I^K + k|)n(|iK-k|)y,±(ri,r2)¥.i(r;,r2), 

^ P K,k 

^±(ri,r2) = ^e'^^ijVi'i?L(r,k), i?i(r, k) = i?z(r, fc)Fi(e,)F;(e,). (80) 

Here the functions G± contain the PDs g±{r) = G±{0, r), and the 1-matrix /(r) = G±{r, oo). 
For 'no polarization' there are the contraction SRs ()79p and /(O) = 1. 



Appendix 3. The case of 'full spin polarization' 

Here C = 1 is considered, i.e. = 0, thus Pj,, 7j,, /j,, = and A''^ = A^. The above 
formulae simplify as 

7i(l|l') = 5<,„.;<5.,t7(ri|r;),7(ri|r;) = p/(rnO, /(O) = 1, (81) 



72(1|1',2|2') = 5+iai\a[,a2\a',)5^,,^5^,,Mr,\r[,r2\r',), Tt^2 = iV(iV - 1), (82) 



p'^ giru) = l{ri\ri,r2\r2), j d^rp[l - g{r)] = 1, 51(00) = 1, 



p^G'(rii/,ri2) 
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7(ri|r'i, r2|r2), ^(0, r^) = gir^), 



j d^ri2 p[f{riv) - G{riy, r^)] = finv), G{riy, 00) = f{riy 
Note giiir) = 0, thus J d^r p[l - gii{r)] = 2/(1 - C) diverges for C ^ 1- 



(83) 
(84) 
(85) 



Appendix 4: The spin structure of 72 
Eq. ()25p agrees term by term with the Eq. (6-13) of Ref. |l^ , where the block-diagonality of 72 
in the representation corresponding to the common eigenstates of the 2-body spin-operators 
Szi + Sz2 and (§1 + §2)^: 

72 = {GiiDiiiG[i + Gi-iDii_iC[_i + GiqDiiqG[q) + GqoDqqqGqq + {Gi()Diq()Gqq + GqoDqioG'iq) 

(86) 



Here Csm, C'gj^j, Dss'm are abbreviations for C^m^ C'^m'; ^S5'A/(ri, r2\r[, r^), respectively. 
For the position matrices it is e.g. Dm = •yii, -Dii-i = 7^, Z^no = 7a-- The triplet /singlet 
spin functions or (Clebsch-Gordan coefficients) 

Cll = ^crit'^<T2T' Cl^l = 6^-^l6cr2l, 

give the corresponding spin matrices of Eqs. ()22j) and according to 

= CuC'ii, 5ii = Ci-iC[_i, (5+ = CioC^o, 

Table I. The 2-matrices of the spin-polarized HEG and their properties (cf. Sec. 5). 
The 6ss and 6^ etc. are defined in Eq.(j221)- The operators D, P, Q, C, and Tr are defined 
after Eq.(j2^. The lines t| and || are for spin-parallel electron pairs, the following lines 
for spin-antiparrallel electron pairs. The 2-matrices depend on ri|r']^, r2|r2, the PDs of 
column D7 depend on = |ri — r2|. In column C7 the 1-matrices f^, 03^4 as well as 
/ = K/t + /i)' / = K/t - fl)^ " = + "4), a = i(a3 - 04) appear, they depend on 
rm = |ri - r;|. Note /^(O) = f\{0) = /(O) = 1 and 03(0) = 04(0) = a{0) = a (which is a 
function of and (), hence /(O) = and (5(0) = 0. 
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Table II. The HF parts of the cumulant expansion (cf. Sec. 8). Note 
af^{r) = a^^{r) = /3^^(r), a^^ = 6^^, a^^ = 0. For 'no interaction' the cumulant 
parts vanish, only the HF parts remain, which simplify furthermore as (3g{r) = jffgi^), 
/3» = §/r(r),thus 6° = f, 6° = §. 
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Table III. The cumulant 2-matrices, being size-extensively normalized (cf. Sec. 8). The 
index 'c' means 'connected' or 'cumulant'. The cumulant normalization TtXss/Ns = Cg 
measures the correlation strength, cf. Eg. ()2()|1 . 
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Table IV. This is Table I, simplified for the case of the spin-unpolarized HEG. In App. 
2 it is shown : 7^ = 7_ = 7^ — 7^, thus also Ag = Qa — Qp, and a = f /2, a = 1/2. Besides, 
the HF parts contain (3^ = Pa = P, thus h^ = ha = h, cf. Eqs.jTHI), 



7 


D7 


P7 


Q7 


C7 


Tr7 


7p 


= (TTTT) 




-7p 


7p 


f/(f-l) 


f(f-l) 


la 


= (TiiT) 


(D'ga 


~l'a 


7a 


p fN 
2J 2 




ia 


= (inT) 




-la 


I'a 


-If 


N 
2 


7+ 


= 7a - 7a 


m9a + Ag) 


+7+ 


7+ 


f/(f +1) 


f(f + 1) 


7- 


= 7a + 7a 


m9a-^9) 


-7- 


7- 


f/(f-l) 


f(f-l) 



[1] K. S. Singwi and M. P. Tosi, in Solid State Physics, ed. by H. Ehrcnrcich, F. Scitz,and D. Turn- 
bull, Academic, New York, 1981, Vol. 36, p. 177; P. Pulde, Electron Correlation in Molecules 
and Solids, 3rd ed.. Springer, Berlin, 1995; N. H. March, Electron Correlation in Molecules 
and Condensed Phases, Plenum, New York, 1996; N. H. March (ed.). Electron Correlation in 
the Solid State, Imperial College Press, London, 1999; W. Kutzelnigg, "Theory of Electron 
Correlation", in J. Rychlewski (ed.), Explicitly Correlated Wave Functions in Chemistry and 
Physics. Theory and Applications, Kluwer, Dordrecht, 2003, p. 3. 

[2] A. Gonis, N. Kioussis, and M. Ciftan (Eds.), Electron Correlations and Materials Properties, 
Kluwer/Plenum, New York, 1999. 

[3] A. Gonis, N. Kioussis, and M. Ciftan (Eds.), Electron Correlations and Materials Properties 
2, Kluwer/Plenum, New York, 2002. 

[4] D.M. Ceperley, Phys. Rev. B 18, 3126 (1978); D.M. Ceperley and B.J. Alder, Phys. Rev. Lett. 
45, 566 (1980). 

[5] J. P. Perdew and Y. Wang, Phys. Rev. B 46, 12947 (1992); 56, 7018 (1997). 

[6] G. Ortiz and P. Ballone, Phys. Rev. B 50, 1391 (1994); 56, 9970 (1997); G. Ortiz, M. Harris, 

and P. Ballone, Phys. Rev. Lett. 82, 5317 (1999). 
[7] P. Gori-Giorgi, F. Sacchetti, and G.B. Bachelet, Phys. Rev. B 61, 7353 (2000); 66, 159901(E) 

(2002). 

[8] P. Gori-Giorgi and J. P. Perdew, Phys. Rev. B 66, 165118 (2002). 

[9] P. Ziesche, phys. stat. sol. (b) 232, 231 (2002); P. Gori-Giorgi and P. Ziesche, Phys. Rev. B 
66, 235116 (2002). 
[10] N. H. March, Phys. Rev. 110, 604 (1958). 

[11] W. Macke and P. Ziesche, Ann. Physik (Leipzig) 13, 25 (1964). 

[12] E. R. Davidson, Reduced Density Matrices in Quantum Chemistry, Academic, New York, 1976. 
[13] J. Cioslowski (ed.), Many-Electron Densities and Reduced Density Matrices, Kluwer/Plenum, 
New York, 2000. 

[14] A. J. Coleman and V. L Yukalov, Reduced Density Matrices. Coulson's Challenge, Springer, 
Berlin, 2000; A. J. Coleman, Phys. Rev. A 66, 022503 (2002); A. Beste, K. Runge, and R. 
Bartlett, Chem. Phys. Lett. 355, 263 (2002). 



[15] E.R. Davidson, Int. J. Quantum Chem. 91, 1 (2003) and refs. therein. 
[16] C. Caccamo, G. Pizzimenti, and M. P. Tosi, Nuovo Cimento 31B, 53 (1976). 
[17] J. C. Kimball, Phys. Rev. A 7, 1648 (1973). 
[18] A. W. Overhauser, Can. J. Phys. 73, 683 (1995). 

[19] P. Gori-Giorgi and J. P. Perdew, Phys. Rev. B 64, 155102 (2001); P. Gori-Giorgi, in Ref. Q, 

p. 379; M. Corona, P. Gori-Giorgi, and J. P. Perdew, unpublished. 
[20] B. Davoudi, M. Polini, R. Asgari, and M.P. Tosi, Phys. Rev. B 66, 075110 (2002); 

|cond-mat /0206456| B. Davoudi, R. Asgari, M. Polini, and M.P. Tosi, unpublished. 
[21] P. Ziesche, Phys. Rev. B 67, 233102 (2003), 

[22] P. Ziesche, K. Pernal, and F. Tasnadi, phys. stat. sol. (b) 239, 185 (2003). 
[23] P. Ziesche and F. Tasnadi, submitted to Ann. Physik (Pobierowo Proceedings, ed. by J. 
Cioslowski) . 

[24] The general parameter theorem dE/dX = (dH/dX) for bound states of H is due to P. 
Giittinger, Z. Phys. 73, 169 (1932), cf. his Eq.(ll). It is later given by W. Pauh, in Handbuch 
der Physik, edited by H. Geiger and K. Scheel (Springer, Berlin, 1933), Band XXIV, Teil 1, 
p. 162; reprinted in Handbuch der Physik, edited by S. Fliigge (Springer, Berlin, 1958), Band 
V, Teil 1, p. 83. A can be any parameter (coupling constant, nuclear coordinate, ...). For a 
coupling constant A and in lowest order of perturbation theory, the theorem was given already 
by E. Schrodinger, Ann. Physik (Leipzig) [4], 80, 437 (1926), cf. his Eq.(7"). The theorem is 
also implicitly contained in Eq.(28) of the paper by M. Born and V. Fock, Z. Phys. 51, 165 
(1928). - Hellmann and later Feynman explicitly referred to the special case of nuclear coordi- 
nates within the Born-Oppenheimer approximation, leading to the 'Hellmann- Feynman' forces 
upon nuclei: H. Hellmann, Einfiihrung in die Quantenchemie, Deuticke, Leipzig, 1937, pp. 61, 
285 (the original russian version is of 23 October 1936: G. Gel'man, Quantenchemie, ONTI, 
Moscow and Leningrad, p. 428); R.P. Feynman, Phys. Rev. 56, 340 (1939). Hellmann stressed, 
that the forces calculated in this way sensitively depend on the electron density used therein. 
For variationally determined densities cf. S.T. Epstein, The Variation Method in Quantum 
Chemistry, Academic Press, New York, 1974. - At degeneracies the theorem is discussed by 
O.E. Alon and L.S. Gederbaum, Phys. Rev. B 68, 033105 (2003). - The Hehmann-Feynman 
theorem has been generalized to Gamow states by P. Ziesche, K. Kunze, B. Milek, J. Phys. A 
20, 2859 (1987). 



[25] M. Isihara, Bull. Chem. Soc. Jpn. 58, 2472 (1985). 

[26] J. Vannimenus and H. F. Budd, Phys. Rev. B15, 5302 (1977). 

[27] Z. Qian and G. Vignale, Phys. Rev. Lett. 88, 056404 (2002); M. Polini and M. P. Tosi, Phys. 

Rev. B 63, 045118 (2001); C. Caccamo, G. Pizzimenti, M. Parinello, and M. P. Tosi, Lett. 

Nuovo Cimento 11, 156 (1974). 
[28] J. C. Kimbah, J. Phys. A 8, 1513 (1975). 
[29] H. Yasuhara and Y. Kawazoe, Physica A 85, 416 (1976). 
[30] P.-O. Lowdin, Adv. Chem. Phys. 2, 207 (1959). 

[31] M. B. Ruskai, J. Math. Phys. 11, 3218 (1970); K. Yasuda, Phys. Rev. A 63, 032517 (2001). 

[32] P. Ziesche, in 0, p. 361; P. Ziesche, J. Tao, M. Seidl, and J. P. Perdew, Int. J. Quantum 
Chem. 77, 819 (2000); P. Ziesche, in Ref. [l^l, p. 33. - For finite systems, particle number 
fluctuations in fragments have been discussed by P. Fulde, cf. Ref. jl|, p. 157, further refs. are 
in Ref. 3- 

[33] A. Rajagopal, J. C. Kimball, and M. Banerjee, Phys. Rev. B 18, 2339 (1978). 
[34] For the cumulant expansion of reduced densities and reduced density matrices cf. e.g. P. 
Ziesche, in jl^, p. 33. 

[35] P. Ziesche, Phys. Lett. A 195, 213 (1994); Int. J. Quantum Chem. 60, 1361 (1996); in Q, p. 

361; M. Levy and P. Ziesche, J. Chem. Phys. 115, 9110 (2001). 
[36] The function p{k) = 3k'^u(k) appears also in the mentioned fluctuation-correlation analysis of 

the spin-unpolarized HEG 22] and in Ref. 37 1. 
[37] P. M. W. Gill, A. M. Lee, N. Nair, and R. D. Adamson, J. Mol. Struct. (Theochem) 506, 303 

(2002). 

[38] P. Ziesche, Int. J. Quantum Chem. 90, 342 (2002); in Ref. 3l> P- 307 

[39] M. Nakata, H. Nakatsuji, M. Ehara, M. Fukuda, K. Nakata, and K. Fujisawa, J. Chem. Phys. 

114, 8282 (2001); cf. also H. Nakatsuji in Ref. |]j], p. 85. 
[40] D. A. Mazziotti, Phys. Rev. A 65, 062511 (2002); cf. also Ref. Q, P- 139. 



Captions of the Figures 



Fig.l. Dimensionless 1-matrix /s(r) of the spin-polarized HEG for different values of the 
spin-polarization ( and the interaction strength r^. r in units of k^^. Upper panel: ideal 
Fermi gas, s =1, and increasing (. Middle panel: ideal Fermi gas, ( = 1/3, and s =t, i- 
Lower panel: ( = 1/3, s =^ , and increasing r^. 

Fig.2. Momentum distributions ns{k,rs,() of the spin-polarized HEG according to 
Eq. ()18j) for the spin-polarization ( = 1/3 and the interaction-strength values = 1,5, 10. 
k in units of kp. 

Fig.3. Kinetic correlation energy tcorr(^s) in a.u. according to Perdew/Wang js^ (full 
line) and from ns{k, r^, () of Eq.lfTHj) via t = +1^ (the dots). 

Fig.4. Non-idempotency matrices l3s{r) and l3a{r) of the spin-polarized HEG according 
to Eq. ()19j) . r in units of kp^. 1st and 2nd panel: ideal Fermi gas, 3rd and 4th panel: 
increasing interaction strength r^. 
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